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Abstract 



We calculate the rates for the charged lepton flavour violating decays 4 — > Ijj, 
t — > in, t — ^ £t], t — y Irj ' , [i~ — > e~e + e~, the six three body leptonic decays 
r~ —7- tjt^i^ and the rate for \i — e conversion in nuclei in the Littlest Higgs model 
with T-Parity (LHT). We also calculate the rates for Kl $ — > ^e, $ — > 7r°/xe and 
Bd,s ~^ £%ky We find that the relative effects of mirror leptons in these transitions 
are by many orders of magnitude larger than analogous mirror quark effects in 
rare K and B decays analyzed recently. In particular, in order to suppress the 
{j, — >■ e7 and yT — > e~e + e~ decay rates and the fj, — e conversion rate below 
the experimental upper bounds, the relevant mixing matrix in the mirror lepton 
sector Vfft must be rather hierarchical, unless the spectrum of mirror leptons is 
quasi-degenerate. We find that the pattern of the LFV branching ratios in the 
LHT model differs significantly from the one encountered in the MSSM, allowing 
in a transparent manner to distinguish these two models with the help of LFV 
processes. We also calculate (g — 2) M and find the new contributions to below 
1 • 10~ 10 and consequently negligible. We compare our results with those present 
in the literature. 



Note added 



An additional contribution to the Z penguin in the Littlest Higgs model with T-parity 
has been pointed out in [1,2], which has been overlooked in the present analysis. This 
contribution leads to the cancellation of the left-over quadratic divergence in the cal- 
culation of some rare decay amplitudes. Instead of presenting separate errata to the 
present work and our papers [3-6] partially affected by this omission, we have presented 
a corrected and updated analysis of flavour changing neutral current processes in the 
Littlest Higgs model with T-parity in [7]. 



1 Introduction 



Little Higgs models [8-10] offer an alternative route to the solution of the little hierarchy 
problem. One of the most attractive models of this class is the Littlest Higgs model 
with T-parity (LHT) [11], where the discrete symmetry forbids tree- level corrections 
to electroweak observables, thus weakening the electroweak precision constraints [12]. 
In this model, the new gauge bosons, fermions and scalars are sufficiently light to be 
discovered at LHC and there is a dark matter candidate [13]. Moreover, the flavour 
structure of the LHT model is richer than the one of the Standard Model (SM), mainly 
due to the presence of three doublets of mirror quarks and three doublets of mirror 
leptons and their weak interactions with the ordinary quarks and leptons. 

Recently, in two detailed analyses, we have investigated in the LHT model AF = 2 
[14] and AF = 1 [3] flavour changing neutral current (FCNC) processes, like particle- 
antiparticle mixings, B — > X s ^, B — > X s £ + £~ and rare K and B decays. The first 
analysis of particle-antiparticle mixing in this model was presented in [15] and the FCNC 
processes in the LH model without T-parity have been presented in [16-20]. 

The most relevant messages of the phenomenological analyses in [3,14,17,20] are: 

• In the LH model without T-parity, which belongs to the class of models with 
constrained minimal flavour violation (CMFV) [21,22], the new physics (NP) effects 
are small as the NP scale / is required to be above 2 — 3 TeV in order to satisfy 
the electroweak precision constraints. 

• In the LHT model, which is not stringently constrained by the latter precision 
tests and contains new flavour and CP- violating interactions, large departures from 
the SM predictions are found, in particular for CP-violating observables that are 
strongly suppressed in the SM. These are first of all the branching ratio for Kl — > 
7r°uu and the CP asymmetry in the B s — > ifxp decay, but also BriK^ — > 
7r°£ + £~) and Br(K + — > tt + uu). Smaller, but still significant, effects have been 
found in rare B s d decays and AM Sjd . 

• The presence of left-over divergences in AF = 1 processes, that signals some 
sensitivity to the ultraviolet (UV) completion of the theory, introduces some the- 
oretical uncertainty in the evaluation of the relevant branching ratios both in the 
LH model [20] and the LHT model [3]. On the other hand, AF = 2 processes and 
the B — > Xg'j decay are free from these divergences. 

Now, it is well known that in the SM the FCNC processes in the lepton sector, like 
U — > £jl and /i~ — > e~e + e~, are very strongly suppressed due to tiny neutrino masses. 
In particular, the branching ratio for /i — > e'-f in the SM amounts to at most 10 -54 , to 
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be compared with the present experimental upper bound, 1.2 • 10~ n [23], and with the 
one that will be available within the next two years, ~ 10~ 13 [24]. Results close to the 
SM predictions are expected within the LH model without T-parity, where the lepton 
sector is identical to the one of the SM and the additional 0(v 2 /f 2 ) corrections have 
only minor impact on this result. Similarly the new effects on (g — 2)^ turn out to be 
small [25]. 

A very different situation is to be expected in the LHT model, where the presence 
of new flavour violating interactions and of mirror leptons with masses of order 1 TeV 
can change the SM expectations up to 45 orders of magnitude, bringing the relevant 
branching ratios for lepton flavour violating (LFV) processes close to the bounds available 
presently or in the near future. Indeed in two recent interesting papers [26,27], it has 
been pointed out that very large enhancements of the branching ratios for — > ij'j and 
t —7- fin are possible within the LHT model. 

The main goal of our paper is a new analysis of li —> ijj, r — > [in, and of other 
LFV processes not considered in [26,27], with the aim to find the pattern of LFV in this 
model and to constrain the mass spectrum of mirror leptons and the new weak mixing 
matrix in the lepton sector Vue, that in addition to three mixing angles contains three 
CP-violating phases 1 . In particular we have calculated the rates for /i~ — > e~e + e~ and 
the six three body leptonic r decays r~ — > ^i^l^, as well as the fi — e conversion rate in 
nuclei. We have also calculated the rates for Kls — > Kl,s ~^ ?rV e an d Bd )S — > ^dj 
that are sensitive to flavour violation both in the mirror quark and mirror lepton sectors. 
Finally we calculated (g — 2) M that has also been considered in [26,27]. 

Our analysis confirms the findings of [26,27] at the qualitative level: the impact of 
mirror leptons on the charged LFV processes — > £j^ and r — > fin can be spectacular 
while the impact on (g — 2)^ is small, although our analytical expressions differ from 
the ones presented in [26,27]. Moreover, our numerical analysis includes also other LFV 
processes, not considered in [26,27], where very large effects turn out to be possible. 

While the fact that in the LHT model several LFV branching ratios can reach their 
present experimental upper bounds is certainly interesting, it depends sensitively on 
the parameters of the model. One of the most important results of the present paper 
is the identification of correlations between various branching ratios that on the one 
hand are less parameter dependent and on the other hand, and more importantly, differ 
significantly from corresponding correlations in the Minimal Supersymmetric Standard 
Model (MSSM) discussed in [29-32]. The origin of this difference is that the dominance 
of the dipole operators in the decays in question present in the MSSM is replaced in the 

^ detailed analysis of the number of phases in the mixing matrices in the LHT model has recently 
been presented in [28]. 
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LHT model by the dominance of Z°-penguin and box diagram contributions with the 
dipole operators playing now a negligible role. As a consequence, LFV processes can 
help to distinguish these two models. 

A detailed analysis of LFV in the LHT model is also motivated by the prospects in the 
measurements of LFV processes with much higher sensitivity than presently available. 
In particular the MEG experiment at PSI [24] should be able to test Br(fi — > e'y) at 
the level of 0(1O~ 13 — 10 -14 ), and the Super Flavour Factory [33] is planned to reach a 
sensitivity for Br(r — > fi'y) of at least (9(10~ 9 ). The planned accuracy of SuperKEKB 
of (9(10~ 8 ) for r — y /ry is also of great interest. Very important will also be an improved 
upper bound on /x — e conversion in Ti. In this context the dedicated J-PARC experiment 
PRISM/PRIME [34] should reach the sensitivity of (9(10~ 18 ), i. e. an improvement by six 
orders of magnitude relative to the present upper bound from SINDRUM II at PSI [35]. 

Our paper is organized as follows. In Section 2 we briefly summarize those ingre- 
dients of the LHT model that are of relevance for our analysis. Section 3 is devoted 
to the decays £j — > £jj with particular attention to fi — > ej, for which a new stringent 
experimental upper bound should be available in the coming years. In Section 4 we 
calculate the branching ratio for r — > fiir and other semi-leptonic r decays for which 
improved upper bounds are available from Belle. In Section 5 we analyze the decays 
/i~ — > e~e + e~, r~ — > fi~fi + fi~ and r~ — > e~e + e~ . In Section 6 we calculate the fi — e 
conversion rate in nuclei, and in Section 7 the decays Kl,s A*e and K^g —> n°fie. 
In Section 8 we give the results for B^ s /^e, re, r/i and in Sections 9 and 10 for 
t~ — > e~fi + e~, /i~e + /i~, /i~e + e~, e~/! + /i~. In Section 11 we calculate (g — 2)^. A de- 
tailed numerical analysis of all these processes is presented in Section 12. In Section 13 
we analyze various correlations between LFV branching ratios and compare them with 
the MSSM results in [29-32]. Finally, in Section 14 we conclude our paper with a list of 
messages from our analysis and with a brief outlook. Few technical details are relegated 
to the appendices. 

2 The LHT Model and its Lepton Sector 

A detailed description of the LHT model and the relevant Feynman rules can be found 
in [3]. Here we just want to state briefly the ingredients needed for the present analysis. 

2.1 Gauge Boson Sector 

The T-even electroweak gauge boson sector consists only of the SM electroweak gauge 
bosons W^, Zl and A^. 



3 



The T-odd gauge boson sector consists of three heavy "partners" of the SM gauge 
bosons 

W§, Z H , A H , (2.1) 
with masses given to lowest order in v/ f by 

M WH =gf, M ZH =gf, M Ah = 9 -^. (2.2) 

All three gauge bosons will be present in our analysis. Note that 

tan#w M Wh 
m a h = —j=-M Wh ~ -jj- , (2.3) 

where 9w is the weak mixing angle. 
2.2 Fermion Sector 

The T-even sector of the LHT model contains just the SM fermions and the heavy top 
partner T+. Due to the smallness of neutrino masses, the T-even contributions to LFV 
processes can be neglected with respect to the T-odd sector. We comment on the issue 
of neutrino masses in the LHT model in Appendix A. 

The T-odd fermion sector [36] consists of three generations of mirror quarks and 
leptons with vectorial couplings under SU(2)l x U(l)y. In this paper, except for Kl,s 
fie, Kl,s — > vr°/ie, Bd, s — > U^j and r — > tix ,tr\,tr{ , only mirror leptons are relevant. We 
will denote them by 

©' 

with their masses satisfying to first order in v/f 

m m = m m , m H2 = m m , m m = m m ■ ( 2 - 5 ) 





1. 1 




0k) 







2.3 Weak Mixing in the Mirror Lepton Sector 

As discussed in detail in [15], one of the important ingredients of the mirror sector is the 
existence of four CKM-like [37] unitary mixing matrices, two for mirror quarks and two 
for mirror leptons: 

V Hu , V Hd , V m , V Hv . (2.6) 

They satisfy 2 

V Hu V H d = Vqkm , V Hv Vm = Vp MNS , (2.7) 

2 Note that it is Vqkm but Vp MNS appearing on the r.h. s., as the PMNS matrix is defined through 
neutrino mixing, while the CKM matrix is defined through mixing in the down-type quark sector. 
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where in Vpmns [38] the Majorana phases are set to zero as no Majorana mass term has 
been introduced for right-handed neutrinos. The mirror mixing matrices in (2.6) param- 
eterize flavour violating interactions between SM fermions and mirror fermions that are 
mediated by the heavy gauge bosons W H , Zh and Ah- The notation in (2.6) indicates 
which of the light fermions of a given electric charge participates in the interaction. 

Thus Vhi, the most important mixing matrix in the present paper, parameterizes the 
interactions of light charged leptons with mirror neutrinos, mediated by W^j, and with 
mirror charged leptons, mediated by Zh and Ah- Feynman rules for these interactions 
can be found in [3]. Vhv parameterizes, on the other hand, the interactions of light 
neutrinos with mirror leptons. 

In the course of our analysis of charged LFV decays it will be useful to introduce the 
following quantities [i = 1,2,3): 



Xi — V H£ v Hl i 



(re) 



V'*2eT tit 
Hi V H£ ? 



v Ht v m i 



(2.8) 



that govern /i — > e, r — > e and r — > \i transitions, respectively. 

We also recall the analogous quantities in the mirror quark sector [i = 1, 2, 3) 



£i - V Hd V Hd > t,i — V Hd V Hd 1 £i — V Hd V Hd 1 l 2 - y J 



.(d) 



,(.s) 



that we will need for the analysis of the decays K L;S — > A^e, K L S — > 7r°/ie and B d s — > £i£j. 

Following [28], we parameterize Vhi i n terms of three mixing angles 8^ and three 
complex phases 5f ■ as a product of three rotations, and introducing a complex phase in 
each of them 3 , thus obtaining 





/l 










J. 

°23 






g 

~ S 23 e 







S 23 e 



'<5.',., 



23 



i 

/ 



13 







1 

\-s[ 3 e iS ™ 



'13 



/ 



/ 



u 12 

s[ 2 e lSl2 



Performing the product one obtains the expression 
/ 



Vi 



Hi 



S 12 C 23 e 12 



\ S 12 S 23 e 



C 12 C 13 
- c\ 2 S 23 s\ 3 e l ^~ S ^ 
J J „iS' 



s\ 2 c\ 3 e ~ia 

S 12 S 23 S 13 e ^ <5l3_<5 * 2 ~' 523 ' > 



s 12 e 



12 



4* o\ 


V 

(2.10) 



i(s 12 +s 23 ) - c\ 2 c 23 s\ 3 e^ 



C 12 C 23 

-4 2 si 3 e i5 ^ - s[ 2 4 3 s[ 3 e l{s ^- s ^ 



s 13 e 

'23 C 13 C 
I I 
C 23 C 13 



(2.11) 

As in the case of the CKM matrix the angles d\ ■ can all be made to lie in the first quadrant 



with < 5{ 2 , 5 23 , 5 l3 < 2rr. The matrix Vhv is then determined through Vh u = VHtVi 



PMNS- 



3 Note that the two additional phases in Vhi have nothing to do with the possible Majorana nature 
of neutrinos. 
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2.4 The Parameters of the LHT Model 



The new parameters of the LHT model, relevant for the present study, are 

/, m e H1 , m e H2 , m e m , 9[ 2 , d{ 3 , 6 e 23 , 5{ 2 5{ 3 5 e 23 , (2.12) 

and the ones in the mirror quark sector that can be probed by FCNC processes in K 
and B meson systems, as discussed in detail in [3,14]. 

The determination of the parameters in (2.12) with the help of LFV processes is 
clearly a formidable task. However, if the new particles present in the LHT model are 
discovered once LHC starts its operation, the parameter / will be determined from M\y H , 
Mz H or Ma h - Similarly the mirror lepton masses m e Hi will be measured. 

The only remaining free parameters among the ones listed in (2.12) will then be 0f- 
and 8fp which can be determined once many LFV processes have been measured. 

3 £j~f in the LHT Model 

3.1 Preliminaries 

In [14] we have shown how one can obtain the branching ratio Br(B — y X^) in the LHT 
model directly from the B — y X s 7 and b — > s gluon decays in the SM by simply changing 
the arguments of the two SM functions D' (x) and E' (x) and adjusting properly various 
overall factors. The explicit formulae for these functions are given in Appendix B. 

Here we will proceed in an analogous way. We will first derive Br(fi — > ej) in the SM 
for arbitrary neutrino masses from the calculation of BriB — y X s '~f) in the same model. 
This will allow us to obtain in a straightforward manner Br(fi — > ej) in the LHT model, 
when also some elements of the Br(B — y X s ~f) calculation in this model in [14] are taken 
into account. The generalization to r — y fi'y and r — y will be automatic. 

The current experimental upper bounds for /x —y ej, r — > /ry and r — y ej are given 
by [23, 39] 4 

Br(n -»• erj) < 1.2 ■ HT 11 , (3.1) 
Br{r ->■ m) < 1.6 • 10 -8 , Br{r -> &y) < 9.4 ■ 10" 8 . (3.2) 

3.2 fi -> ej in the SM 

The diagrams for the B — y Xg'y and fi — y e-y decays in the SM are shown in Figs. 1 
and 2, respectively. In gauge also the corresponding diagrams with Goldstone bosons 

4 The bounds in [39] have been obtained by combining Belle [40] and BaBar [41] results. 
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Figure 1: Diagrams contributing to B —y X s ~f in the SM. 




Figure 2: Diagrams contributing to /i — >• e7 m i/ie 5M. 

contribute. The diagram with the photon coupled directly to the internal fermion line, 
present in B — > X s j, is absent in the case of fi — y due to the neutrino charge neutrality. 

In the case of the B — y decay the function D' (x) resulting from the diagrams in 
Fig. 1 can be decomposed as follows 

D' (x) = TO*)] Abclian + [^(x)] tripIe , (3.3) 

where the first term on the r. h. s. represents the sum of the first and the last diagram 
in Fig. 1 and the second term the second diagram. 

The inspection of an explicit calculation of the b — y 57 transition in the 't Hooft- 
Feynman gauge gives 

PoOOlAWtan = Wd - Qu)E' Q (x) , (3.4) 

with Qd and Q u being the electric charges of external and internal fermions, respectively. 
Setting Qd = —1/3 and Q u = 2/3 in (3.4) and using (3.3) we find 

[£>o(z)Wie = D' (x) + ~E' (x) . (3.5) 

This contribution is independent of fermion electric charges and can be directly used in 
the fj, — y &y decay. 

Turning now our attention to the latter decay we find first from (3.4) 

[«L-2^)> (3-6) 

as Q u = and Qd = — 1 in this case. The final result for the relevant short distance 
function in the case of \i — y is then given by the sum of (3.5) and (3.6). Denoting 
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this function by H(x\) we find 



H(x l v ) = D' (xl) - -K(x 



2 



3" 0V vJ ' ~ u \M W/ 
The generalization of the known SM result to arbitrary neutrino masses reads then 

3a 



x, 



m. 



(3.7) 



Br(fi — > cy) 



SM 



2tt 



.1 , 



(3.{ 



with Vjj being the elements of the PMNS matrix. 
Now, in the limit of small neutrino masses, 



(3.9) 



and we confirm the known result 



Br(/i ->■ e7) S M 



3a 
32^ 



El/ .V*-x i 
'er {ii^v 



(3.10) 



Assuming that the angle #13 of the PMNS matrix is not very small so that Am 



atm 



dominates the branching ratio, we find 

Br([i -> e7) SM ^ 0.015 a 
where Am^ tm < 3.2 • 10" 3 eV 2 and s l3 < 0.2 [42] 



Am 2 , 



M 2 W 



< 10 



-51 



(3.11) 



3.3 n -> e~f in the LHT Model 

The diagrams contributing to \x — > &y in the LHT model are shown in Fig. 3. We 
show only the contributions from the mirror fermion sector as the T-even sector gives a 
negligible contribution. Note that the heavy scalar triplet $ does not contribute at this 
order mv/f (see [3,14] for details). 

Let us write the resulting branching ratio as follows 



3a 
2^ 



Br(fx ->■ e7) L HT = tti \^w H + ^z H + A Ah \ , 



(3.12) 



with the different terms representing the W^, Zjj and An contributions. 
Defining 



m 



Vi 



Hi 



y[ = a i/i with a 



tan 9w 



16.6, 



(3.13) 
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Figure 3: Diagrams contributing to /x — > cy in the LET model. 



and including the factor 



we find, using (3.8) 



M w H 
lv 2 



1«2 

4/5' 



(3.14) 



(3.15) 



with H defined in (3.7) and xf in (2-8). 

The neutral gauge boson contributions can directly be deduced from (4.10) of [14]. 
Including a factor 3 that takes into account the difference between the electric charges 
of quarks and leptons, we obtain from the last two terms in (4.10) of [14] 



A 



lv 



772 x 



(pe) 



4/ 
lv 2 



E 



(tie) 



\nivi 



(3.16) 
(3.17) 



Finally, adding the three contributions in (3.15)-(3.17), we find using (3.12) 



Br{n -»■ e7) LHT 



3a 
2^ 



D 



I pe 1 2 
odd ' 



with 



pe 

-^odd ^ j2 



E 



X 



(pe) 



D' ( yi ) - 7 -E' (y t ) - ^E'M) 



(3.18) 



(3.19) 



with xt defined in (2.8), yi in (3.13) and D' Q , E' given in Appendix B. The formulae 
(3.18) and (3.19) represent the main result of this section. 
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Let us next compare our result with the one of Goyal [26], which has been obtained 
without including the interferences between A Wh , A Zh and Aa h , whose omission cannot 
be justified. His results for |Avy ff | 2 , |A^| 2 and (A^j 2 agree with ours, provided the 
latter contribution is divided by 25 in (15) of [26], once one assumes the straightforward 
relation Fa h (Za h ) = F ZH (aZ Zll ) between the A H and Z H short distance functions. Note 
that the short distance functions Fj(zj) in [26] still contain mass independent terms that 
disappear in the final expressions after the unitarity of the Vhi matrix has been used. For 
this reason we have removed such terms from the functions D' (x) and E' (x). Finally, in 
contrast to [26] , we find that the contributions of the scalar triplet $ are of higher order 
in v/f. 

In a recent paper by Choudhury et al. [27] a new analysis of £j — > Iff has been 
presented in which the interference terms have been taken into account. Unfortunately, 
the formulae in that paper are quite complicated thus preventing us from making an 
analytic comparison. On the other hand, we have performed a numerical comparison, 
finding significant differences between our results and those in [27]. 

3.4 t — y e~f and r — ► [ij in the LHT Model 

The branching ratios for r — > e^ and r — > /17 can easily be found in analogy to fi — > e^y. 
We find 

Br{r ery) = ^Br(r~ v T e~v e ) \D' ^ , (3.20) 

Br(r ^) = ^Br(r~ v rV Tv^ \D' Z\ 2 , (3.21) 

where D'J dd and D'J dd can be obtained from (3.19) by replacing (fie) with (re) and (rfi), 
respectively. Furthermore [43] 

Br(r~ -> v T e~v e ) = (17.84±0.05)% , Br(r~ -> i/ T //~P M ) = (17.36±0.05)% . (3.22) 

4 Semi-leptonic r Decays 

Recently the Belle [44] and BaBar [45] collaborations presented improved upper bounds 
for the decays r — > IP (P = n, rj, rj'), which have been combined [39] to 

Br(r -> fin) < 5.8 ■ 10~ 8 , Br(r -> firf) < 5.1 ■ 10" 8 , Br(r -> fir]') < 5.3 • 10~ 8 , 

(4.1) 

Br(r -> en) < 4.4 • 10" 8 , Br(r ->■ erf) < 4.5 • 10~ 8 , Br(r -> erf') < 9.0 ■ 10" 8 , 

(4.2) 
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W H 



W h 



Wh' 




• n 



Wh\ 



■ Wh Zh , Ah ' 



>Zh,Ah Zh,Ai 



>Zh, Ah 



Figure 4: Diagrams contributing tor—)- fin in the LHT model. Similar diagrams, but 
with d quarks in the final state, also contribute. 



thus increasing the interest in investigating these branching ratios in the LHT model. 

The branching ratios for these semi-leptonic decays have already been estimated 
within the LHT model by Goyal [26], but with the rough approximation consisting in the 
use of the effective Hamiltonian for K° — K° mixing. Here we will study the semileptonic 
decays in question with a more accurate approach, having at hand the recent analysis of 
rare K and B decays in the LHT model [3]. 

The diagrams for r — > fin are shown in Fig. 4. As tt° has the following flavour 
structure 

r0 _uu- dd ^ ^ 



7T ' 



V2 



there are two sets of diagrams, with uu and dd in the final state. The corresponding 
effective Hamiltonians can directly be obtained from [3]. They involve the short- distance 
functions X^ A and — Y^£j for uu and dd, respectively. Taking into account the opposite 
sign that is conventionally chosen to define the two short distance functions, the effective 
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Hamiltonian that includes both sets of diagrams is simply given as follows 



H 



{ x Za( uu )v-A - Y^ d (dd) v „ A ) (fir) 



V-A 



(4.4) 



sin 2 6 W 

Here X^ dd and Y^ dd have the same structure as the functions calculated in [3] in the 
context of rare K and B decays. Adapting them to the lepton sector we find: 

(J ua (y 2 , z) - J u ^y u z)) + (J ua (ys, z) - J uU ( Vl , z)) 
(J dd (y2, z) - j*{ Vu z)) + x ^ (J d3 (ys, z) - J dd ( yi , z)) 



^odd 



1 odd 



(4.5) 
(4.6) 



where 



64? 



S , 



dd 



+4 (G( Vi , z; Z H ) + GM, z'; A H ) + G 2 (y h z; rjfj 
y t S odd + F dd ( yi , Z] W H ) 

A(G( yu z; Z h ) + Gite, z'; A H ) - G 2 ( Vll z; 77)) 
(4tt/) 2 



(4.7) 



1 . yf 



log 



(4.8) 
(4.9) 



with the functions F ua , F dd , G, Gi and G 2 given in Appendix B, the leptonic variables 
fji and y[ defined in (3.13) and the analogous variables for degenerate mirror quarks 5 
given by 



m 



g 2 
H 



1 

z = a z , T] = — 
a 



As 



(0\(uu) v _ A \n°) = -(0\(dd) v _ A \n°) 



V2 



where F n ~ 131 MeV is the pion decay constant, we find 



Br(r — > /x7r) 



G 2 nrF i m s r 

U F U r n'"'T l T | yT/J, i y-TjU|2 
^odd 1 odd I ' 



(4.10) 
(4.11) 

(4.12) 



128vr 3 sin 4 9 W 

with r T and m T being the lifetime and mass of the decaying r, and neglecting suppressed 
pion and muon mass contributions of the order 0(m 2 T /m 2 .) and (9(m 2 /m 2 ). The branch- 
ing ratio for the r — > en decay can be obtained very easily from (4.12) by simply replacing 



J The limit of degenerate mirror quarks represents here a good approximation, as the box contributions 
vanish in the limit of degenerate mirror leptons and, consequently, the inclusion of mass splittings in 
the quark sector represents a higher order effect. 
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(r/i) with (re). The generalization of (4.12) to the decays r — » [ir\ and r — > /jt/ is quite 
straightforward too, although slightly complicated by mixing in the 77 — 7/ system. 

The understanding of 77 — 77' mixing has largely improved in the last decade mainly 
thanks to the formulation of a new mixing scheme [46,47], where not one but two angles 
are introduced to relate the physical states (77, 77') to the octet and singlet states (77s, 770), 
as 

\r]) = cos 9 8 \r] 8 ) - sin 9 \r] ) , \r]') = sm9 8 \r] 8 ) + cos9 \r] ) , (4.13) 

with 



\Vs) = 77g(M + \dd) ~ 2\ss)) 



\ r]0 ) = ^=(\uu) + \dd) + \ss)). (4.14) 



In this mixing scheme, four independent decay constants are involved. Each of the 



two physical mesons (P 
components, defined by 



77,77'), in fact, has both octet (a = 8) and singlet (a = 0) 



(0|(?yg)v-A|P(p)> = ^ 



(4.15) 



where the SU(3) generators A a satisfy the normalization convention Tr[A a A 6 ] = 25 ab . 
They are conveniently parameterized [46] in terms of the two mixing angles (9 8 , 9 ) and 
two basic decay constants (F 8 ,F ), as 



778 770 

7] 7] 
t r ] ' t r ] ' 



F 8 cos9 8 -F sm9 
F 8 sin 9 8 F cos 9q 



(4.16) 



Working in this mixing scheme and generalizing the expression for the r — > fiir 
branching ratio in (4.12), one obtains 

2 

Fn - 

7(T 



Br(r -)■ /J77) - 



128vr 3 sin 4 # 



ir 



cos 9 8 F 8( - 



odd J 



gSin^ (J 



odd 



9 ' 
z 2 odd; 



/J77 



128tt 3 sin 4 W 



sin 0s F> 



x/3 ^ 



8 / yTjl 



odd 



odd 



9 ' 
z 2 odd; 



(4.17) 



We conclude this section noting that the function S^dd in (4.9) represents a left-over 
singularity that signals some sensitivity of the final results to the UV completion of the 
theory. This issue is known from the study of electroweak precision constraints [11] 
and has been discussed in detail in the context of the LH model without T-parity [20] 
and, more recently, also with T-parity [3]. We refer to the latter paper for a detailed 
discussion. Here we just mention that, in estimating the contribution of these logarithmic 
singularities as in (4.9), we have assumed the UV completion of the theory not to have a 
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complicated flavour pattern or at least that it has no impact below the cut-off. Clearly, 
this additional assumption lowers the predictive power of the theory. In spite of that, we 
believe that the general picture of LFV processes presented here is only insignificantly 
shadowed by this general property of non-linear sigma models. 



5 \i — > e e + e , r — > fi and r — > e e + e 

Next, we will consider the decay fi~ — > e~e + e~, for which the experimental upper bound 
reads [48] 

Br{n~ -> e"e + e-) < 1.0- 1(T 12 . (5.1) 

This decay is governed, analogously to the h — > sfi + fi~ transition, analyzed in the LHT 
model already in [3], by contributions from 7- and Z°-penguins and by box diagrams. 
However, the fact that now in the final state two identical particles are present does 
not allow to use directly the known final expressions for b — > sfi + fi~, although some 
intermediate results from the latter decay turned out to be useful here. Also the general 
result for fi~ — > e~e + e~ obtained in [49], which has been corrected in [29,31], turned out 
to be very helpful. 

Performing the calculation in the unitary gauge, where we found the contribution 
from the Z°-penguin to vanish [3], we find for the relevant amplitudes from photon 
penguins and box diagrams 6 : 

G e 2 1 

A' = yf^^odd [e(pi)(m^CT Q /3^(l + 7 5 )V(p)] ® [e(p2h a e(p 3 )} 

-(Pi****), (5-2) 

G F e 2 



r2 -^odd [e(Pi)7a(l - 7s)Mp)] ® [e{p2)l a e{p 3 )} - (pi ^ p 2 
Gp a 



W2STX 2 



(5.3) 



A ox = 2^ , q o e dd [e( Pl )7a(l - 75MP)] O [e(p 2 )7 a (l - iMPs)} ■ (5-4) 

y 2 2TT Sin U\\r 

The function -D odd is given in (3.19), while the functions Y^ dd and Z^ dA can easily be 
obtained from those calculated in [3]. The analogy with the b — > sfi + fi~ decay, together 
with the observation that the fi" — > e~e + e~ decay in question involves only leptons in 



'Following [49], our sign conventions are chosen such that T-L e s is determined from —A. 
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both the initial and final states, allow us to write' 



e,odd 



Era 2 U^y^-J^Vi) 



i=l 



+ x? e) E i^i 2 Jdd ^ v*) - Jdd ^ v* 



(5.5) 



with J given in (4.8). Following a similar reasoning we can write for the Z^ dd function 



^odd 



X^ £) {Z odd (y 2 ) - Z odA (yi)) + xT ] { Z odA{yz) - Z odd {yx)) 



where 



Zodd(yi) = C dd(yi) + ^-Dodd(z/i) • 



(5.6) 



(5.7) 



The explicit expressions of the C odd and D odd functions are given in Appendix B 8 . 
Here, we just note that as a consequence of the charge difference between the leptons 
involved in f/T —> e~e + e~ and the quarks involved in b — > sfi + /i~, D odd in (5.7) differs 
from the analogous function found in [3]. 

Comparing these expressions to the general expressions for the amplitudes given 
in [31,49], we easily obtain r(/i~ — > e~e + e~). Normalizing by r(/i" — > e~v e v IJ ), we find 
the branching ratio for the decay fi~ — > e~e + e~ to be 

r(/i~ — > e~e + e~) 



Br(n — > e e e~ 



T(/i ->■ e u e u^ 



Q- 



71" 



6 Z odd 



+ 3 Re (Z^DZIY) + \D' T d \ Z ( log ^ - ^ 



11 



2 sin 4 6^' e ' odd l sin 2 ^ 



sin 2 6 



For t — > /i we make the following replacements in (5.2)-(5.8): 



->■ m T , 



(5.8) 



(5.9) 



so that, in particular, Y^ dd is now present. Furthermore, in (5.8) the normalization 
T(/i~ — > e~u^u e ) is replaced by r(r~ — > ^rv T v^), so that the final result for Br{r~ — > 



The subscript e of Y^ dd denotes which of the SM charged leptons appears on the flavour conserving 
side of the relevant box diagrams. 

8 Note that the functions C dd and -Dodd arc gauge dependent and have been calculated in the 't Hooft- 
Feynman gauge. However, the function Z^ dd is gauge independent, so that it can be used also in the 
unitary gauge calculation above. 
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fi fi + fi ) contains an additional factor Brir —t \i z^-z^). In the case of r — > e e + e 
the replacements in (5.5)-(5.8) amount only to 

(fie) — > (re) , m M — >• m r , (5.10) 

having now Y™ dd , and in (5.8) r(//~ — >■ e'v^v^) is replaced by r(r _ — > e~v T u e ) and 
an additional factor Br(r~ —> e~u T u e ) appears. In doing this we neglect in all three 
expressions m eyfl with respect to m T . 



6 /x — e Conversion in Nuclei 

Similarly to the decays \i — > and fi~ — > e~e + e~ , stringent experimental upper bounds 
on /j, — e conversion in nuclei exist. In particular, the experimental upper bound on \i — e 
conversion in ^Ti reads [35] 

R(fiTi -> eTi) < 4.3 • 10^ 12 , (6.1) 

and the dedicated J- PARC experiment PRISM/PRIME should reach a sensitivity of 
C(10~ 18 ) [34]. 

A very detailed calculation of the \x — e conversion rate in various nuclei has been 
performed in [50], using the methods developed by Czarnecki et al. [51]. It has been 
emphasized in [50] that the atomic number dependence of the conversion rate can be 
used to distinguish between different theoretical models of LFV. Useful general formulae 
can also be found in [49]. 

We have calculated the fi — e conversion rate in nuclei in the LHT model using the 
general model-independent formulae of both [49] and [50] . We have checked numerically 
that, for relatively light nuclei such as Ti, both results agree within 10%. Therefore, we 
will give the result for /i — e conversion in nuclei derived from the general expression 
given in [49], as it has a more transparent structure than the one of [50]. 

Following a similar reasoning as in the previous section, we find from (58) of [49] 

r(/iX^eX) = ga 5 ^|F(g)| 2 mJ 



y-fie y/ze 

Z (4^ d e d + D' Z) - (2Z + N)—lf- + (Z + 2A0- 



sin 9 W sin 9 W 



(6.2) 



where X£ dd and Y^ dd are obtained from (4.5) and (4.6) by making the replacement 
(r/i) — >■ (fie), and D' odd and Z^ dd are given in (3.19) and (5.6), respectively. Z and iV 
denote the proton and neutron number of the nucleus. Z e g has been determined in [52] 
and F(q 2 ) is the nucleon form factor. For X = ^Ti, Z e g = 17.6 and F(q 2 ~ —mf,) — 0.54 
[53]. 
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The fi — e conversion rate i?(/xX — > eX) is then given by 

R(P X -> eX) = ^ ^ , (6.3) 

^ capture 

with r^ pture being the /i capture rate of the element X. The experimental value is given 
by rapture = (2-590 ± 0.012) • 10 6 s" 1 [54]. 

In our numerical analysis of Section 12 we will restrict ourselves to fi — e conversion 
in ^Ti, for which the most stringent experimental upper bound exists and where the 
approximations entering (6.2) work very well. For details, we refer the reader to [49,50, 
53]. 

7 K L S fie and K L S -¥ 7r°/xe in the LHT Model 

The rare decay Kl — > fie is well known from the studies of the Pati-Salam (PS) 
model [55], where it proceeds through a tree level leptoquark exchange in the t-channel. 
The stringent upper bound on its rate [56] 

Br(K L -»■ yLe) = Br(K L /i+e") + Br(K L yTe + ) < 4.7 ■ 10" 12 , (7.1) 

implies the mass of the leptoquark gauge boson to be above 10 3 TeV [57]. By increasing 
the weak gauge group of the PS model in the context of the so-called "Petite Unifica- 
tion" [58] and placing ordinary fermions in multiplets with heavy new fermions, and not 
with the ordinary SM fermions as in the PS model, it is possible to avoid tree level con- 
tributions to Kl — y fie so that the process is dominated by the box diagrams with new 
heavy gauge bosons, heavy quarks and heavy leptons exchanged. The relevant masses of 
these particles can then be decreased to (9(1 TeV) without violating the bound in (7.1). 

In the SM the decay Kl /-te is forbidden at the tree level but can proceed through 
box diagrams in the case of non-degenerate neutrino masses. Similarly to fi — > e7 it is 
too small to be measured. 

In the LHT model as well, Kl — > fie appears first at one loop level. It proceeds 
through the diagrams shown in Fig. 5 that should be compared with very similar dia- 
grams in Fig. 1 of [14] contributing to particle-antiparticle mixing in the T-odd sector 
of the LHT model. The main difference is the appearance of leptons in the lower part of 
the box diagrams instead of quarks, leading to a different structure of the loop functions 
once the unitarity of the matrices Vhcl and Vhi has been used. 

The effective Hamiltonian for Kl — > /ie can be directly obtained from (3.11) of [14] 
by removing the QCD factor 772, appropriately changing the CKM-like factors and the 
arguments of the short distance functions and multiplying by 2 to correct for a different 
combinatorial factor. 
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>Wh Z h ,Aj 



>Zh,Ah Zh,Ai 



>Zh,Ah 



f fl H 



Figure 5: Diagrams contributing to Kl — > fi + e in the LHT model. Similar diagrams 
describe K L — > /i~e + . 



Taking into account that both the /i + e and /i e + final states are experimentally 
detected, the relevant effective Hamiltonian reads 



+ X ( r ) \sd) V -A^e) v _ A +h.c, (7.2) 



where x\ an d Hi are defined in (2.8) and (3.13), and ^■ JX) and Z{ belong to the quark 
sector and are defined as 



JK) 



— v Hd V Hd > 

The short distance functions read 



Hi 



(7.3) 



F H (zi, yj) = F(zi, yj] W H ) + G(zi, yj\ Z H ) + Ai(zi, yy, Z H ) + A 2 (z h yy Z H ) , (7.4) 

where the different contributions correspond to "WW", "ZZ", U AA" and U ZA" dia- 
grams, respectively. Explicit expressions for the functions F, G, Ai and A 2 can be found 
in Appendix B. 

Using the unitarity of the Vna and Vm matrices we find 



G% 2 v 2 



xtki K) R(z 2 , z u yi , y 2 ) + x rct'R(^ *i, vi, vs) 



+ X2^3 R i z 3, Z U Vi, V2) + xT'iT > R{z2, zt, y u y 3 ) {sd) v _ A (efj 1 ) V - A 



+ 



X ( r r ^ K) R(z 2 , Zl , Vl , y 2 ) + x f*> &>R(z 3 , z u Vl , y 3 ) 



+ X { 2 £) ^ 3 K) R(z 3 ,z 1 ,y 1 ,y 2 ) + xT' £2^(^2,^1,2/1,2/3) (sd) V - A (fie) 



IV— A 



+ h.c. 



(7.5) 
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where 



R(zi,Zj,y k ,yi) = F H (Zi,yi) + F H (zj,y k ) - F H (zj,yi) - F H (zi,y k ) . 



(7.6) 



We remark that the relevant operators differ from (c?7 M e)(/27 M s) present in the PS 
model which is characteristic for transitions mediated by leptoquark exchanges. There- 
fore, the branching ratio for Kl — > /ie in the LHT model is most straightforwardly 
obtained from the one for Kl — V /i + /i~ in the SM, after the difference between /i + /i~ 
and /ie has been taken into account. 

With the help of (XI.44) and (XXV. 1) of [59] we easily find 



Br(K L /ie) = Br(K L /i + e~) + Br(K L n~e 



(7.7) 



128vr 4 " Wl / 4 



Br(K + /iV 



r{K L 



t(K+) |K. 



i,j=2,3 



(7.8) 



Here [43,60] 

Sr(A'+ -> /i+z/) = (63.44 ± 0.14)% 



t(K+) 



4.117±0.019 



IK 



0.225 ±0.001 . 

(7.9) 

The corresponding expression for Br(Ks —> /ie) is obtained from (7.7) through the 
replacements [43] 



t{K l ) k t(K s 



t(K+) t(K 



(7.229 ±0.014) ■ 10" 



(7.10) 



Due to t(Ks) <C t(Kl), the branching ratio Br(K$ — > /ie) is expected to be typically 
by two orders of magnitude smaller than Bt(Kl — > /ie) unless Ini(£ t - ) ^> Re(£^). 

The decay Kl — > ir°fie in the LHT model is again governed by the effective Hamil- 
tonian in (7.5). This time it is useful to perform the calculation of the branching ratio 
in analogy with Kl — > 7r uu [61]. Removing the overall factor 3 in BriK^ ^vv) 
corresponding to three neutrino flavours, we find 



Br{K L -> ttV) = Br(K L -> 7r°/i + e") + Br{K L -> 7r°fi~e + ) 

128^ /4^ K ^vr/i ^ r( ^ +)|Ks|2 



M=2,3 



(7.11) 
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Here [43], 

Br(K + -> ir°fi + v) = (3.32 ± 0.06)% . (7.12) 

We note that this time lm(^ K ^) instead of Ke(^ K ^) enters, coming from the difference 
in sign between the relations 

(7r°\(ds) v . A \K°) = -(n°\(sd) v . A \K°) (7.13) 

and 

(0\(ds) V - A \K°) = +(0\(sd) v _ A \K°) . (7.14) 

The corresponding expression for Br(Ks —> ir fie) is obtained from (7.11) through 
the replacements 

T M^ T M' Im «, m )^Re«r>). (T.15) 

We would like to emphasize that in deriving the formulae for Bt{Kl,s fi-e) and 
BriyKi s ~~ >" vr°/ie) we have neglected the contributions from CP violation in K — K° 
mixing (indirect CP violation). For instance in the case of BriK^ — > ir°fie) only CP 
violation in the amplitude (direct CP violation) has been taken into account. The indirect 
CP violation alone gives the contribution 

Br(K L -» 7rV) ind ' = Br(K s ->■ ir°fie)\e K \ 2 (7.16) 

and needs only to be taken into account, together with the interference with the con- 
tribution from direct CP violation, for Im^- ) <C Re^-^). The latter case, however, 
is uninteresting since it corresponds to an unmeasurably small branching ratio. This 
should be contrasted with the case of Kl — > ir°e + e~, where the indirect CP violation 
turns out to be dominant [62]. The origin of the difference is that photon penguins, 
absent in Kl — > 7r°/ie, are present in Kl — > ir°e + e~ and the structure of "H e ff is rather 
different from (7.5). Moreover, while the estimate of indirect CP violation to Kl — > ir°fie 
in the LHT model can be done perturbatively, this is not the case for Kl — > 7r°e + e~, 
where the matrix elements of the usual AF = 1 four quark operators have to be taken 
into account together with renormalization group effects at scales below M\y. 

In summary the branching ratios for Bt(Kl,s — > fi-e) and Bt(Kl,s — ► ^°fie) in the 
LHT model can be calculated fully in perturbation theory and are thus as theoretically 
clean as Kl — > tt uu. As seen in the formulae (7.7), (7.10), (7.11) and (7.15) above, 
Bt{Kl — > fie) and Br(Ks — > k° fie) are governed by Re(£p^), while Br {Kg — > fie) 
and Bt(Kl — > it fie) by lm(^ K ^). Neglecting the indirect CP violation, we find then 
that in the so-called "i^-scenario" of [3], the T-odd contributions to BriK^ — > fie) and 
Br(Ks —> n° fie) are highly suppressed, while Br(Ks —> fie) and Bt{Kl —> n°fie) are 
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generally non-vanishing. The opposite is true in the case of the "S s -scenario" in which 
only Br(Ki — > fie) and Br(K$ — > 7T°fie) differ significantly from zero. 

This discussion shows that the measurements of Br(Ki — > fie) and Br(Ki — > n fie) 
will transparently shed some light on the complex phases present in the mirror quark 
sector, and from the point of view of the LHT model, the measurement of BriK^ — > 
Tr°fie) at the level of 10~ 15 will be a clear signal of new CP-violating phases at work. 



8 B d , s ->> /xe, B d , s ->> re and B d:S — >► r/x 



The decays of neutral 5-mesons to two different charged leptons proceed similarly to 
the Kl —> fie decay discussed in Section 7. The effective Hamiltonian describing these 
processes receive again contributions only from box diagrams. For the Bd — > fie decay, 
it reads 



n eS {B d fie) 



r 2 

u f /If 2 

32tt 2 Wl * 2 



7 2 



Y^^FH^y,) \x { r\bd)v-A{efi)v- 



+ X [ r ] \bd) V - A {fie)v- 



■1) 



with t^f* = VuV^Hd- Using the unitarity of the Vud and Vm matrices, it becomes 
n cS (B d ^fie) = Yl tFRfatuVuVi) \xt e) (pd)v-A(efi) V - A 



f 

+ X i r > (bd) v . A (fie) v . 



(8.2) 



with R(zi, Zj,i/k,yi) being the combination of short distance functions defined in (7.6). 

The effective Hamiltonians describing the remaining decays have a similar structure 
and can easily be derived from "H c ff (-E^ — > fie) in (8.2) through the following replacements 





—> fie) : 


d d) -+ } , 


H- c R{B d 


-> re) : 


xT } x { ; e) > 




-> re) : 


A d ) . A s ) O e ) . Are) 

Si ' Si ) Xj * Xj 


"Hcfi {B d 


->■ Tfi) : 


x { r ] x [ r ] , 


'Hes(B s 


-»• rfi) : 


A d ) _> As) A^ e ) (t/j.) 
Si ' Si ) Xj * Xj 



(8.3) 



The effective Hamiltonians for the B^ s decays are simply given by the hermitian conju- 
gates of the corresponding expressions in (8.2) and (8.3). 
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In calculating the corresponding branching ratios, it is important to observe that 
while in Kl — > fie the decaying meson, is a mixture of flavour 

eigenstates, the decaying B^ iS mesons are instead flavour eigenstates (e.g. B d = bd). For 
this reason, in Br(Ki — > fie) the two conjugate contributions of % e ff combine together, 
while in B ds decays only one contribution enters, with the conjugate one describing B ds 
decays. With this consideration in mind, starting from Br[K B — > fie) in (7.7), one can 
straightforwardly find the following expressions 



Br(B d -> fie) 



Br(B d ->■ fie) = Br(B d -> /i + e~) + Br(B d -> /i~e~ 



512vr 4 |K b | 2 / 4 r(^ 



U) 



M=2,3 



»J=2,3 



Br(B s ->■ /ie) 



Br(B B ->> /xe) = 5r(5 s -»■ ,u + e~) + £r(£ s -»■ /i"e H 
G\M^ l v*r(B s )M Bs Fl 



5127i*\V ub \ 2 f*T(B+)M Bd F* d 



5.5) 



M=2,3 



i,j=2,3 



£r(5 d -> re) 



£r(£ d -> re) = Br(B d -> r+e") + Br(S d -> r"e + ) 
G 2 F M^ ^ r (S d ) 



512vr 4 |K b | 2 / 4 r(^ 



-£r(£ + -»■ r + z/ r ) 



(8.6) 



M=2,3 



£ Q xi R(zi,zi,yuVj] 

i,j=2,3 



Sr(5 s ->■ re) = £r(£ s ->> re) = Br(5 s -> r + e~) + Br{B s -»■ r~e + ) 



5127r 4 |Kb| 2 / 4 t(B + ) M Bd F B 



■Br(B + -> r + z/ T 



5.7) 



£ ^ i *i R ( z i> z l>yi>Vj) 



i,J=2,3 
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Br(B d -> Tfi) 



Br(B d -»• r/i) = 5r(5 rf -> r+/0 + 5r(B d -> rV 
^ v'r(B d ) 



5127r 4 |K b | 2 / 4 ^(5+ 



-£r(£ + -»■ r + z/ r ) 



(8.8) 



i,i=2,3 



£r(£ s ->> r/x) = 5r(5 s -> r/x) = Br(B s -»■ r + /0 + £r(£ s r~/i" 



512vr 4 |K b | 2 / 4 ^-(5+)M Sd F j 



(8.9) 



Analogously to Br(K^ — y fie), we have chosen to normalize the branching ratios in 
question introducing the branching ratio of a B + leptonic decay. This normalization will 
become helpful once Br(B + — y ^v^) and Br(B + — y t + u t ) are experimentally measured 
with sufficient accuracy. Currently the measurements for Br(B + — y t + v t ) read 



Br(B + -> r + i/ r ) 



(1-7918^) • 10- 
(O.SSj™^) ■ 10' 



[63] 
[64] 



5.10) 



while the experimental upper bound on Br(B + — y fi + Vfj) is given by [65] 

Br(B + ->■ // + ^) < 6.6 ■ 10" 6 (90% C.L.) . 



(8.11) 



Therefore, in our numerical analysis in Section 12 we will use the central values of the 
theoretical predictions for these decays, given by 



Br(B+ -> n+u„ 



G 



V ub \ 2 F 2 Bd M Bd mlr(B + ) = (3.8 ± 1.1) ■ 10" 7 



Br(B+ -> r + u T ) = -^\V ub \ 2 F 2 B M Bd m 2 T T{B + ) = (1.1 ± 0.3) ■ 10' 



57T 



(8.12) 
(8.13) 



where the relevant input parameters are collected in Table 1 of Section 12. 



9 t — Y e /x + e and r — > /x e + /x 

These two decays are of AL = 2 type and are very strongly suppressed in the SM. In 
the LHT model they proceed through the box diagrams in Fig. 6. 
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Figure 6: Diagrams contributing to t — > e jj, + e in the LHT model. 



The effective Hamiltonians for these decays can be obtained from AB = 2 processes, 
that is from (3.11) of [14]. In the case of r~ — > e~/i + e~ we find 



U 



g 2 v 2 

u f m 2 — 



XiXj F H(yi, %)(er)y-A(e/i)v_A 



(9.1) 



i,3 



with the function Fjj given in (7.4). The additional factor 4 relative to (3.11) of [14] 
results from the different flavour structure of the operator involved and the two identical 
particles in the final state. 

The corresponding effective Hamiltonian for r~ — > fi~e + fi~ is obtained by simply 
exchanging e and /i, with xf^ — xf^ ■ 

The relevant branching ratios can be found by comparing these two decays to the 
tree level decay — > v T e~v ei for which the effective Hamiltonian reads 

G p 



UcE = -^y T r) v _ A {ev e )y_ A 



and yields the decay rate 

T(r~ — > v T e~v e 
For r~ — > e~fi + e~ we find then 



G\m\ 
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Br(r — > e fi + e 
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1927T 3 V 167T 2 J f 



h3 



(9.2) 



(9.3) 



(9.4) 



where we neglected m Me with respect to m T . We have included the factor 1/2 to take 
into account the presence of two identical fermions in the final state. 

The branching ratio for t~ — > ^Te + yT is obtained from (9.4) by interchanging fi -h- e. 



10 t — y \x e + e and r —ye /x + /x 

These decays have two types of contributions. First of all they proceed as in t~ — >■ 
fi~fi + fi~ and t~ — > e~e + e~ through AL = 1 penguin and box diagrams. As this time 
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there are no identical particles in the final state, the effective Hamiltonians for these 
contributions can be directly obtained from the decay B — > X s £ + £~. Let us derive 
explicitly the effective Hamiltonian for r _ — > fi~e + e~. The generalization to t~ — > 
e~ fi + fi~ will then be automatic. 

As the QCD corrections are not involved now, only three operators originating in 
magnetic photon penguins, Z°-penguins, standard photon penguins and the relevant 
box diagrams have to be considered. Keeping the notation from B — > X s fi + fi~ but 
translating the quark flavours into lepton flavours these operators are 

Q 7 = -^m r /ia Q/3 (l + l5 )rF af5 , (10.1) 

that enters, of course with different external states, also the fi — > decay, and 

Q 9 = (fLT) V - A (ee) v , Qio = (flr)v-A(ee) A • (10.2) 

The effective Hamiltonian is then given by 

U cS (t- /i-e+e") = -5| [C/Q 7 + C/Q 9 + C^Qio] . (10.3) 

The Wilson coefficient for the operator Q 7 can easily be found from Section 3 of the 
present paper and Section 7 of [3]. We find 

C t / = -\d' Z, (10.4) 

with D'JJ^ obtained from (3.19) by replacing (fie) with (rfi). 

The Wilson coefficients of the operators Q$ and Qiq receive not only contributions 
from AL = 1 7-penguin, Z°-penguin and box diagrams, but also from AL = 2 box 
diagrams, as in Section 9. For C^ 1 and C[q we can then write 

cr = £c7 , c% = %fis , (10.5) 

CI/ = -r^r - 4^dd - A r , , = — ^f- + A T , , (10.6) 

Sin U\y sm "w 

with the functions Yj^ dd and Z T odA obtained from (5.5) and (5.6) by replacing (fie) by 
(rfi). A Tfl represents the additional AL = 2 contribution which is not present in the 
case of b — > s£ + £~ and will be explained below. As there are no light fermions in the 
T-odd sector, the mass independent term present in Cg in the case of b — > s£ + £~ in (X.5) 
of [59] is absent here. Effectively this corresponds to setting 77 = 1 in the latter equation 
and of course removing QCD corrections. 

As already mentioned, also the AL = 2 diagrams shown in Fig. 7 contribute to this 
decay. The corresponding effective Hamiltonian can be obtained from (9.1) through 
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Figure 7: Diagrams of AL = 2 type contributing to r —*■(/, e + e in the LHT model. 



the obvious replacements of local operators, removing the symmetry factor 2 and the 
following change in the mixing factors: 



Aj A.j 



(-) (Me)* 
Aj Aj 1 



so that we find 



A 
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16 sin 2 9 W f 2 



i.j 



(10.7) 



(10.8) 



Effectively the presence of the diagrams in Fig. 7 introduces corrections to the Wilson 
coefficients Cg and C w in (10.6). As the relevant operator has the structure (V — A) <g) 
(V — A), the shifts in Cg and Cio are equal up to an overall sign. 

Finally, introducing 



{Pe+ +Pe~f 



R T »(s) 



|r(r- -4 fi-e+e- 
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r(T ^z/ r ) 



(10.9) 



and neglecting m e with respect to m T we find for the differential decay rate R T,1 (s) 

,2 



^(S) = — (1-5) 



1 + 2s) ( |C 9 r/ f + l^' 2 
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+4(1 + ;) |C 7 T/ f + 12 Re ( C^CZ") 



;io.io) 



The branching ratio is then given as follows: 



5r(r -»■ /i e + e ) — Br(r -> /x / R TfJ, (s)ds. 



(10.11) 



The branching ratio for r — > e can easily be obtained from the above expres- 

sions by interchanging \i +>• e, where = • 
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For quasi-degenerate mirror leptons the AL = 1 part clearly dominates as the GIM- 
like suppression acts only on one mirror lepton propagator, whereas it acts twice in the 
AL = 2 case. Moreover, in the latter case the effective Hamiltonian is quartic in the Vm 
couplings, whereas it is to a very good approximation quadratic in the case of AL = 1. 
As these factors are all smaller than 1, quite generally AL = 2 contributions will then 
be additionally suppressed by the mixing matrix elements. Consequently, the AL = 1 
part is expected to dominate and the shift A rAt can be neglected. On the other hand, for 
very special structures of the Vm matrix, the double GIM suppression of AL = 2 with 
respect to AL = 1 contributions could be compensated by the Vm factors. Therefore it 
is safer to use the more general expressions given above. 

11 (g - 2) M in the LHT Model 

The anomalous magnetic moment of the muon = (g — 2)^/2 provides an excellent test 
for physics beyond the SM and has been measured very precisely at the E821 experiment 
[66] in Brookhaven. The latest result of the (g — 2) Collaboration of E821 reads 

a* xp = (11659208.0 ± 6.3) ■ 10" 10 , (11.1) 

whereas the SM prediction is given by [67] 

fl SM = a QED + fl cw + a had = ( 11659180 .4 ± 5.1) ■ 10" 10 . (11.2) 

While the QED and electroweak contributions to a^ M are known very precisely [68,69], 
the theoretical uncertainty is dominated by the hadronic vacuum polarization and light- 
by-light contributions. These contributions have been evaluated in [67,70-72]. 

The anomalous magnetic moment a M can be extracted from the photon-muon vertex 
function Y^(p',p) 

u(p')T»(p',p)u(p) = u(p') [^Fviq 2 ) + {v + v'YF M {q 2 )\ u(p) , (11.3) 

where the anomalous magnetic moment of the muon can be read off as 

a M = -2mF M (0). (11.4) 

The diagrams which yield new contributions to a M in the LHT model are shown in 
Fig. 8. They either have a heavy neutral gauge boson [Zn or Ah) and two heavy charged 
leptons £ H (i — 1, 2, 3) or two heavy charged gauge bosons (W|) and one heavy neutrino 
v l H {% = 1, 2, 3) running in the loop. 
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Figure 8: Diagrams contributing to (g — 2) M in the LHT model. 



Calculating the diagrams in Fig. 8 and using the Feynman rules given in [3], the 
contributions of the new particles for each generation i — 1, 2, 3 are found to be: 
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11.6) 



11.7) 



11.8] 



The parameter m Hi in (11.5) and (11.6) denotes the mass of the mirror leptons while Mx 
is the mass of the heavy gauge bosons. We expanded our results in the small parameter 
m^/Mx- Our results in (11.5) and (11.6) for the muon anomalous magnetic moment are 
confirmed by the formulae in [73] for general couplings. 

Replacing the parameters r and b by the more convenient parameter defined in 
(3.13), leads us to the following expressions 
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(11.9) 
(11.10) 
(11.11) 
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m e -- 


= 0.5110 MeV 




r(B d )/r(B + ) = 0.934(7) 




= lUo.oo MeV 




nB s ) 


= 1.466(59) ps 


m T 


= 1.777(J(oj Gev 




t(B+) 


= 1.638(11) ps 


r T = 


: 290.6(10) -10 ps 




M Bd = 


5.2794(5) GeV 


M w 


nn /lor /on\ f *■ T 7" 

= 80.425(38) GeV 




71 JT 

M Bs = 


5.3675(18) GeV [43] 


a = 


1 /I Q*7 
l/iO( 




\v ub \ = 


3.68(14) ■ 10" 3 [75] 


Gp 


= 1.16637(1) ■ 10~ 5 


GeV- 2 




= 1.28 (ChPT) 


sin 2 


6 W = 0.23122(15) 


[43] 




= 1.18(4) 


F Bd 


= 189(27) MeV 




08 = ~ 


22.2(18)° 


Fb s 


= 230(30) MeV 


[74] 


e = - 


8.7(21)° [76] 



Table 1: Values of the experimental and theoretical quantities used as input parameters. 



where the functions Li and L 2 are given in Appendix B. 
Our final result for in the LHT model therefore is 
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While we disagree with [26], we confirm the result of [27] except that according to 
us the factors (V^^V^^i and (V^^V^^ in equations (3.22)-(3.24) of that paper 
should be replaced by |V^| 2 . 



12 Numerical Analysis 
12.1 Preliminaries 

In contrast to rare meson decays, the number of flavour violating decays in the lepton 
sector, for which useful constraints exist, is rather limited. Basically only the constraints 
on Br(fi — > e^y), Br(fi~ — > e~e + e~), i?(/iTi — > eTi) and Br(Ki — > fie) can be mentioned 
here. The situation may change significantly in the coming years and the next decade 
through the experiments briefly discussed in the introduction. 

In this section we want to analyze numerically various branching ratios that we have 
calculated in Sections 3-11. In Section 13 we will extend our numerical analysis by 
studying various ratios of branching ratios and comparing them with those found in the 
MSSM. Our purpose is not to present a very detailed numerical analysis of all decays, 
but rather to concentrate on the most interesting ones from the present perspective and 
indicate rough upper bounds on all calculated branching ratios within the LHT model. 
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To this end we will first set / = 1 TeV and consider three benchmark scenarios for the 
remaining LHT parameters in (2.12), as discussed below. 

In Table 1 we collect the values of the input parameters that enter our numerical 
analysis. In order to simplify the analysis, we will set all input parameters to their 
central values. As all parameters, except for the decay constants Fs d s and the 77 — 77' 
mixing angles, are known with quite high precision, including the error ranges in the 
analysis would amount only to percent effects in the observables considered, which is 
clearly beyond the scope of our analysis. 



12.2 Benchmark Scenarios 

We will consider the following three scenarios: 

Scenario A (red): 

In this scenario we will choose 

Vm = V PMNS , (12.1) 

so that Vhv = 1) an d mirror leptons have no impact on flavour violating observables 
in the neutrino sector, such as neutrino oscillations. In particular we set the PMNS 
parameters to [42] 

sin# 12 = V0.300, sinfl 13 = V0.030 , sin0 23 = ^, 5 13 = 65° , (12.2) 

V2 

which is consistent with the experimental constraints on the PMNS matrix [43]. As no 
constraints on the PMNS phases exist, we have taken <5 13 to be equal to the CKM phase 
and set the two Majorana phases to zero. 

Furthermore, we take the mirror lepton masses to lie in the range 

300GeV<m^ < 1.5 TeV, (i = 1, 2, 3) . (12.3) 



Scenario B (green): 

Here, we take 

V m = V C KM, (12.4) 

so that [75] 

4 = 13°, 4 = 0.25°, 4 = 2.4°, (12.5) 

4 = 0, 4 = 65°, 4 = 0, (12.6) 

and the mirror lepton masses in the range (12.3). 
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Figure 9: Correlation between fi — > cy and fi~ — > e~e + e~ in the scenarios of Section 12.2. 
In the right plot of Scenario C we show the contributions to /i~ — > e~e + e~ from DJ^ 
(purple, lowermost), Z^ dd (orange, middle) andY^ e odd (light-blue, uppermost) separately. 
The shaded area represents the experimental constraints. 

Scenario C (blue): 

Here we perform a general scan over the whole parameter space, with the only restriction 
being the range (12.3) for mirror lepton masses. 

At a certain stage we will investigate the dependence on mass splittings in the mirror 
lepton spectrum. 

In the case of Kl — > fie, B^ s — > fie and similar decays of Sections 7 and 8, the 
parameters of the mirror quark sector enter and the constraints from K and B physics, 
analyzed in [3, 14], have to be taken into account. 

12.3 /x — V e7, fi~ — >- e~e + e~ and fj, — e Conversion 

In Fig. 9 we show the correlation between fi — > cy and /i~ — > e~e + e~ in the scenarios in 
question together with the experimental bounds on these decays. We observe: 
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• In Scenario A the great majority of points is outside the allowed range, implying 
that the Vm matrix must be much more hierarchical than Vpmns in order to satisfy 
the present upper bounds on /i — > and ji~ — > e~e + e~ . 

• Also in Scenario B most of the points violate the current experimental bounds, 
although Vckm is much more hierarchical than Vpmns- The reason is that the 
CKM hierarchy S13 <C S23 <C S12 implies very small effects in transitions between 
the third and the first generation, like r — > e^y, while allowing relatively large effects 
in the [i — > e transitions. Thus in order to satisfy the experimental constraints on 
/1 — > e7 and fi~ — > e~e + e~ a very different hierarchy of the Vm matrix is required, 
unless the mirror lepton masses are quasi-degenerate. 

• In Scenario C there are more possibilities, but also here a strong correlation between 
H — > e7 and \i~ —> e~e + e~ is observed. This is easy to understand, as both decays 
probe dominantly the combinations of Vm elements Xi ■ 

• For Scenario C, we also show the contributions to fi~ — > e~e + e~ from D' odd , Z^ dd 
and Y^ dd separately. We observe that the dominant contributions come from the 
functions Z^ dd and above all Y^ dd , while the contribution of the operator Q 7 , given 
by D' odd , is by roughly two orders of magnitude smaller and thus fully negligible. 
This should be contrasted with the case of the MSSM where the dipole operator is 
dominant. We will return to the consequences of this finding in the next section. 

• We emphasize that the strong correlation between [i — > e^y and fi~ — > e~e + e~ in 
the LHT model is not a common feature of all extensions of the SM, in which 
the structure of pT — > e~e + e~ is generally much more complicated than in the 
LHT model. It is clear from Fig. 9 that an improved upper bound on fi — > &y by 
MEG in 2007 and in particular its discovery will provide important information on 
pT — > e~e + e~ within the model in question. 

Next, in Fig. 10 we show the correlation between the /i — e conversion rate in 22 Ti 
and Br(fi — > ej), after imposing the existing constraints on fi — > &y and fi~ — > e~e + e~ . 
We observe that this correlation is much weaker than the one between /1 — > &y and 
fi~ — > e~e + e~ . Furthermore, we find that the fi — e conversion rate in Ti is likely to be 
found close to the current experimental upper bound, and that in some regions of the 
parameter space the latter bound is even the most constraining one. 

We also show in Fig. 11 the mass splitting Am 12 = \m e H2 — m e Hl \ as a function of m H1 , 
after imposing the existing constraints on fi — > and /i~ — > e~e + e~. We observe that 
for the scenarios A and B considered here, the first two generations of mirror leptons 
have to be quasi-degenerate in order to satisfy the experimental constraints. While in 
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R (/jTi^eTi ) 




Br (M->ey) 
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Figure 10: \x — e conversion rate in ^Ti as a function of Br (fx — > ej), after imposing 
the existing constraints on \i — >• ej and \i~ — > e~e + e~. The shaded area represents the 
current experimental upper bound on R(fiTi — > eTi). 




Figure 11: Mirror lepton mass splitting Ami2 = |^H2 — m Hi\ as a function ofm l H1 , after 
imposing the existing constraints on fi — >• ej and fi~ — > e~e + e~ . The dark points violate 
the constraint on ji — e conversion in Ti, while the light points fulfil also this constraint. 
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Br (r-^jUy) 




Figure 12: Correlation between Br(r — > cy) and Br[r — > fi'y). 

the case of Scenario A, Am 12 ^ 40 GeV is required, Am 12 ^ 60 GeV is sufficient to 
fulfil the constraints in Scenario B, provided the mirror lepton masses are relatively 
small. Generally we find that the allowed mass splittings are larger for smaller values 
of the mirror lepton masses. This is due to the fact that the functions J uu and J dd , as 
defined in (4.7) and (4.8), are monotonously increasing, making thus a stronger GIM 
suppression necessary in the case of large masses. Finally we have studied the impact of 
the experimental bound on R(fiTi — > eTi) on the allowed mass splittings. We observe 
that in particular for large values of the mirror lepton masses, fi — e conversion turns 
out to be even more constraining than fi — > ey and fi~ — > e~e + e~ . 

VIA r — y fi'j and r — )• e'j 

In Fig. 12 we show the correlation between Br(r — > /ry) and Br{r — > cy) in the scenarios 
considered, imposing the experimental bounds on fi — > ey and fi~ — > e~e + e~. We 
observe that they both can be individually as high as ~ 8 • 10~ 10 , but the highest values 
of Brir — > fi^y) correspond generally to much lower values of Brir — \ cy) and vice versa. 
Still simultaneous values of both branching ratios as high as 2 • 10~ 10 are possible. 

12.5 r — > fin, fir], fir]' and r — > /i.7 

In Fig. 13 we show Brir — > fin) and Br(r — > firj) as functions of Br(r — > fi'y), imposing 
the constraints from fi — > &y and /i~ — > e~e + e~. We find that Brir — > fin) can reach 
values as high as 2 • 10~ 9 and Brir — > fit]) can be as large as 7 • 10~ 10 , which is still by 
more than one order of magnitude below the recent bounds from Belle and BaBar. We 
do not show Br{r — > fir/') as it is very similar to Br{r — > /irj). 
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Figure 13: Brir — > fin) (dark-blue) and Br(r — > firj) (light-blue) as functions of Brir — > 
/i 7 ). 

Completely analogous correlations can be found also for the corresponding decays 
r — > e7T, er], erf and r — > e7- Indeed, this symmetry between r — > fi and r — > e systems 
turns out to be a general feature of the LHT model, that can be found in all decays 
considered in the present paper. We will return to this issue in Section 13. 

An immediate consequence of these correlations is that, as in the case of r — >■ fi'y and 
r — > cy, the highest values for r — > fin are possible if r — > en is relatively small, and 
vice versa. Still the corresponding branching ratios can be simultaneously enhanced to 
3 • 10~ 10 . Analogous statements apply tor-y /i(e)?7 and r — > fi(e)r}'. 

12.6 K L — >- fie and K L — >- 7r%te 

In Fig. 14 we show the dependence of Br(Ki — > fie) on Br(fi — > e^y). To this end we 
have chosen one of the sets of mirror quark parameters for which the most spectacular 
effects both in and the K — > nvv decays have been found [3]. We observe that for 
the parameters used here, BriK^ — > fie) is still by two orders of magnitude below the 
current experimental upper bound. However, one can see in Table 2 that Br(Ki — > fie) 
could also be found only one order of magnitude below the current bound. This means 
that large effects in the K — > nvv decays do not necessarily imply also large effects in 
Kl — > fie. 

The effects in Kl — > n°fie are even by roughly two orders of magnitude smaller, as 
can also be seen in Table 2, so that, from the point of view of the LHT model, this decay 
will not be observed in the foreseeable future. 
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1 * ■ — ■ — Br(^ey) 

l.xicr 15 l.xicr 14 l.xio -13 l.xicr 12 l.xicr 11 

Figure 14: Correlation between Br {fx — > e^y) and BriK^ — > fie) for fixed mirror 
quark parameters {m q H1 ~ 332 GeV, m q H2 ~ 739 GeV, m q m ~ 819 GeV, 6f 2 — 95°,#f 3 ~ 
229°, <9| 3 ~ 205°, 5f 2 = 0, 8f 3 ~ 239°, 5 2 d 3 = 0). 

12.7 Upper Bounds 

In Table 2 we show the present LHT upper bounds on all branching ratios considered 
in the present paper, together with the corresponding experimental bounds. In order to 
see the strong dependence on the scale /, we give these bounds both for / = 1000 GeV 
and / = 500 GeV with the range (12.3) for the mirror lepton masses in both cases. We 
observe that the upper bounds on r decays, except for t~ — > fj,~e + /j,~ and r _ — > e~ fi + e~ 
increase by almost two orders of magnitude, when lowering the scale / down to 500 GeV, 
so that these decays could be found close to their current experimental upper bounds. In 
particular, the most recent upper bounds on r — > fin, en [39] could be violated by roughly 
a factor 5, as can be seen from Table 2 of the first version of the present paper. Therefore, 
in deriving the LHT upper bounds for / = 500 GeV, we have taken into account also 
the latter bounds. On the other hand, the bounds on r~ — \ fi~e + fi~ and r _ — > e~ fi + e~ 
are quite independent of the value of /. This striking difference is due to the fact that 
the present lepton constraints are only effective for \i — y e transitions. We also note 
that the upper bounds on some K and B^ s decays, namely Kl — > fie, — > n°fxe and 
Bd,s A*e, result to be lower when the NP scale is decreased to / = 500 GeV. The 
origin of this behaviour is that lowering / the strongest constraints, mainly on K and 
B systems, start to exclude some range of parameters and, consequently, to forbid very 
large values for the branching ratios in question. 

We have also investigated the effect of imposing in addition the constraint R{fiT\ — > 
eTi) < 5 ■ 10~ 12 , which we choose slightly above the experimental value 4.3 • 10~ 12 in 
order to account for the theoretical uncertainties involved. We find that all upper bounds 
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Table 2: Upper bounds on LFV decay branching ratios in the LHT model, for two different 
values of the scale f , after imposing the constraints on fi —> ej and fi~ —> e~e + e~ . The 
numbers given in brackets are obtained after imposing the additional constraint R(fiTi — y 
eTi) < 5 • 10" 12 . For f = 500 GeV, also the bounds on r — > fin, en have been included. 
The current experimental upper bounds are also given. 
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collected in Table 2 depend only weakly on that constraint. This finding justifies that 
we did not take into account this bound in our numerical analysis so far, as it has only 
a minor impact on the observables discussed. 

We would like to stress that the bounds in Table 2 should only be considered as rough 
upper bounds. They have been obtained from scattering over the allowed parameter 
space of the model. In particular, no confidence level can be assigned to them. The 
same applies to the ranges given in Table 3 for the LHT model. 

12.8 (g - 2) M 

Finally, we have analyzed (g — 2) M in the LHT model. Even for the scale / being as low 
as 500 GeV, we find 

aJ HT <1.2-10- 10 , (12.7) 

to be compared with the experimental value in (11.1). We observe that the effect of 
mirror fermions is by roughly a factor of 5 below the current experimental uncertainty, 
implying that the possible discrepancy between the SM value and the data cannot be 
solved in the model considered here. 

13 Patterns of Correlations and Comparison with 
Supersymmetry 

13.1 Preliminaries 

We have seen in the previous section that the branching ratios for several charged LFV 
processes could reach within the LHT model the level accessible to experiments per- 
formed in this decade. However, in view of many parameters involved, it is desirable 
to look for certain correlations between various branching ratios that are less parameter 
dependent than individual branching ratios, and whose pattern could provide a clear 
signature of the LHT model. 

In the case of CMFV in the quark sector useful correlations have been summarized 
at length in [22,83]. In the case of LFV, in a very interesting paper [29], Ellis et 
al. noticed a number of correlations characteristic for the MSSM, in the absence of 
significant Higgs contributions. These correlations have also been analyzed recently 
in [30-32]. In particular, in [30,32] modifications of them in the presence of significant 
Higgs contributions have been pointed out. 

The main goal of this section is a brief review of the correlations discussed in [29-32] 
and the comparison of MSSM results with the ones of the LHT model. We will see that 
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indeed the correlations in question could allow for a transparent distinction between 
MSSM and LHT, which is difficult in high energy collider processes [84]. 



13.2 Correlations in the MSSM 
13.2.1 Dipole Operator Dominance 

In the absence of significant Higgs boson contributions, the LFV processes considered in 
our paper are dominated in the MSSM by the dipole operator with very small contribu- 
tions from box and Z-penguin diagrams. In this case one finds the approximate general 
formulae [29-32] 



Br(£- ^ £j£f£j) a / m\ 



Br(£i £j~f) ' 3tt \ & m 



log Hf-2.7 , (13.1) 



Br & ^ £-£+£-) a I mj. 

' log^- 2.7 . 13.2 

mj. 



Br(£i -> ' 3tt 

Consequently one finds 

£>r(/i~ — > e~e + e~) a 



log-£-2.7) (13.3) 



Br(n — > cy) 3ir \ ml J 162 



Br(r — > cy) Br{r — > fi^) 3ir \ ml J 95 

Bt{t- -)• Br(r~ —> e~fi + fi~) a ( m 2 T \ 1 

log— 5- - 2.7 ~ — , (13.5) 



Br(r — > /jL'y) Br(r — >■ e7) 37r \ y 438 

and [32] 

Br(r~ —7- uTe + e~) Br(r~ — > e~e + e~) 

t - v-4.6. (13.6) 

Br(r~ -» /i~/i + /i") Br(r~ — > e~/i + /i~) 

Moreover, keeping only the dipole operator contribution in i?(/iTi — >■ eTi), we find 

i^Ti ->■ eTi) 



-Br(/i — » e7) 

One also has 



0.7a. (13.7) 



l%^ e *y) <0{a) ( P = ™1, (13.8) 

where the absence of dipole operator contributions to Br(r — > £P) makes the ratios in 
(13.8) significantly smaller than the ones in (13.3)-(13.5). 
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13.2.2 Including Higgs Contributions: Decoupling Limit 



It should be emphasized that the Higgs contributions become competitive with the gauge 
mediated ones, once the Higgs masses are roughly by one order of magnitude smaller 
than the sfermion masses and tan/3 is O(40 — 50). There is a rich literature on Higgs 
contributions to LFV within supersymmetry. One of the earlier references is the one 
by Babu and Kolda [85]. Here we concentrate on the modifications of the correlations 
discussed above in the presence of significant Higgs contributions as analyzed by Paradisi 
[32], where further references can be found (see in particular [30]). 

In the limit of Higgs decoupling, which is achieved in the MSSM, some of the results 
in (13.1)— (13.8) are modified. In particular 



Br(r~ -)■ fi~fi + fi~) 2 Br(r" -)■ e~fi + fi~) 1 



Br(r 9 ' Br(r -»■ erf) 12 

can be much larger than in the case of dipole operator dominance, and 



13.9) 



Br(r — > a e + e ) . Br(r — > e e + e ) . . 
7TT , \ Z0.05, — -j — 4>0.13 13.10 

can be much smaller. Moreover one finds 

Brir — > irj) , Br(r — > fit]) Br(r —> erf) 

1 — 4 5 — 12 

Br(r — > £7) ~ Br(r~ — > n~n + n~) Br(r~ — > e~ji + ji~) 

(13.11) 



13.3 Correlations in the LHT Model 

The pattern of correlations in the LHT model differs significantly from the MSSM one 
presented above. This is due to the fact that the dipole contributions to the decays 
l~ —> Pj^^j an d £~ —> Pjf't^k can De ^lly neglected in comparison with Z°-penguin 
and box diagram contributions. This is dominantly due to the fact that the neutral 
gauge boson (Zjj, Ah) contributions interfere destructively with the contributions 
to the dipole operator functions 5 ^ d , but constructively in the case of the functions 
^fc odd t na t summarize the Z°-penguin and box contributions in a gauge independent 
manner. Moreover, the large tan/3 enhancement of dipole operators characteristic for 
the MSSM is absent here. In this context it is useful to define the ratios 

2 

(13.12) 

that are strongly enhanced in the case of the LHT model for almost the entire space 
of parameters, as seen in Fig. 15 for the case of T Me . Consequently, the logarithmic 



T 



1 j,odd 



-^odd 



10 




Figure 15: T Me as a function of \D'J^ d \ for f = 1 TeV. 

enhancement of dipole contributions seen in (5.8), (13.1) and (13.2) is eliminated by the 
strong suppression of |-D d d | 2 with respect to l^oddl 2 - Similarly, C^, that is governed 
by D'J dd , can be neglected in (10.10), resulting in the absence of large logarithms in the 
decays r~ — > /i~e + e~ and r~ — > e~/i + /i~. 

We note that the ratios Tij depend on the approximation made for the left-over 
singularity, as Y^ odd suffers from this divergence while D*i A does not. However, we 
find that dropping completely the term proportional to S odd) which is certainly not a 
good approximation, amounts to at most 50% changes in the allowed ranges for T^- 
and consequently for the correlations discussed below. Therefore we think that these 
correlations are only insignificantly affected by this general feature of non-linear sigma 
models. 

These findings make clear that the correlations between various branching ratios in 
the LHT model, being unaffected by the logarithms present in (13.1)— (13.5), have a very 
different pattern from those found in the MSSM. 

It turns out that within an accuracy of 3% one can set D' Q z dd to zero in all decays 
with three leptons in the final state. Similarly, the contribution of C$ in r~ — > fi~e + e~ 
and r~ — > e~fi + fi~ can be neglected. Neglecting finally in (10.6), which is good to 
within 20%, allows us to derive simple expressions for various ratios of branching ratios 
that can directly be compared with those listed in (13.1)— (13.6) . 

To this end, we define 



i.odd 



\ri3 
1 fc.odd 



j,odd 



13.13) 
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We find then 



Qj = 3 sin 4 9w \ a>ij | 2 + ~ ~~ 2 sin 2 6*wRe(ajj) . (13.14) 



£?r(/i — > e7) 37r sin 4 

with analogous expressions for the respective r — >■ e and r — )■ fi transitions. 
We also find 

Br(r~ —> e~fi + fi~) a 1 



fir(r e7) " 127rsin 4 #n/ 
Br(r~ —> ji~e + e~) a 1 



fir(r^/i 7 ) " 127rsin 4 ^ 
Finally, (13.15)-(13.17) imply 

Br[j~ — > e~e + e~) c 



T re K e . (13.16) 



T Tfl b*, (13.17) 



-re 



Br(r~ — > e~/i+yu~) 6re 



Br{r~ — >• /i~e+e~) 6 e 



13.18) 
13.19) 



In what follows we restrict ourselves to / = 1 TeV for simplicity. We note that 
although the numerical values given below depend slightly on the size of the scale /, the 
qualitative picture found and discussed remains true independently of that value. 

The ranges for the ratios in question found in the LHT model are compared in Table 
3 with the corresponding values in the MSSM, both in the case of dipole dominance and 
when Higgs contributions are significant. 

While the results in Table 3 are self-explanatory, let us emphasize four striking dif- 
ferences between the LHT and MSSM results in the case of small Higgs contributions: 

• The ratio (13.15) and the similar ratios for r — > e and r — > \i transitions are 0(1) 
in the LHT model as opposed to O(o) in the MSSM. 

• Also the fi — e conversion rate in nuclei, normalized to Br(fi — > ej), can be signifi- 
cantly enhanced in the LHT model, with respect to the MSSM without significant 
Higgs contributions. However the distinction in this case is not as clear as in the 
case of Br(£~ — > iji^ij)/ Br(£i — > £jj) due to a destructive interference of the 
different contributions to /i — e conversion (see (6.2)). On the other hand, in the 
case of the MSSM with significant Higgs contributions, i?(//Ti — > eTi) /Br(fi — > &y) 
is typically much larger than a, so that a distinction from the LHT model would 
be difficult in that case. 
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Table 3: Comparison of various ratios of branching ratios in the LHT model and in the 
MSSM without and with significant Higgs contributions. 

• The "inverted" pattern of the ratios in (13.6) is absent in the LHT model as 
the ratios (13.18) and (13.19) are comparable in magnitude. Moreover they are 
generally very different from the MSSM values. 

• The last finding also implies that while the ratios (13.4) and (13.5) differ roughly 
by a factor of 5 in the case of the MSSM, they are comparable in the LHT model, 
as seen in (13.16) and (13.17). 

In order to exhibit these different patterns in a transparent manner, let us define the 
three ratios 

Br(r~ — > e~e + e~) Br(r~ — > fi~e + e~) 



R2 
R 3 



Br(r~ — > \i~ fJ, + fj,~) Br(r~ — > e~yU + yU") 

Br(r~ — > e~e + e~) Br(r — > /ry) 
Br(r~ — > /i-yU + /i~) Br(r — > ej) 

Br{r~ — » e~fi + fi~) Br(r — > yiry) 



(13.20) 
(13.21) 
(13.22) 



Br{r~ — > /i~e + e - ) Br(r — > cy) 

Note that in the case of a fi -H- e symmetry, these three ratios should be equal to unity. 
This symmetry is clearly very strongly broken in the MSSM, due to the sensitivity of 
the ratios in (13.1) and (13.2) to m e and m M , where one finds 



i?i ~ 20 



R-2 



R* ~ 0.2 



(MSSM) . 



^13.23) 
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On the other hand, in the case of the LHT model the absence of large logarithms allows 
to satisfy the /x -B- e symmetry in question within 30%, so that we find 



0.8 < R x < 1.3, 0.8 < R 2 < 1.2 , 0.8 < R 3 < 1.2 (LHT). (13.24) 

The comparison of (13.23) with (13.24) offers a very clear distinction between these two 
models. 

In the presence of significant Higgs contributions the distinction between MSSM 
and LHT is less pronounced in r decays with /x in the final state. However even in 
this case both models can be distinguished, as seen in Table 3. In particular, the four 
ratios involving Br(£i — > are still significantly smaller in the MSSM than in the LHT 
model, and all decays with electrons in the final state offer excellent means to distinguish 
these two models. 

For the decays r — > £P with P = tt, rj, rj' we find the ranges 

. Br(r £ir) . . Br(r ->> £rj) . . Br(r -)> in') , 

1 < 1 L < 5 5 4 < - — < 2 3 < - — < 2 8 

£?r(r £7) fir(r ->■ £7) £?r(r ->■ £7) 

(13.25) 

As seen in (13.8), the ratios obtained in the LHT model are much larger than in su- 
persymmetry without significant Higgs contributions. On the other hand, if the Higgs 
contributions are dominant, the distinction through the first inequality of (13.11) will 
be difficult. However, we also find 

, Br(r — > ari) , , Br(r — > en) 
0.7 < — \ ^'l <1.3, 1.1 < — 1 ^-<l-8, 13.26 

that could be distinguished from the corresponding results in (13.11). 

In summary we have demonstrated that the LHT model can be very transparently 
distinguished from the MSSM with the help of LFV processes, while such a distinction 
is non-trivial in the case of high energy processes [84]. We consider this result as one of 
the most interesting ones of our paper. 



14 Conclusions 

In the present paper we have extended our analysis of flavour and CP- violating processes 
in the LHT model [3, 14] to the lepton sector. 

In contrast to rare K and B decays, where the SM contributions play an important 
and often the dominant role in the LHT model, the smallness of ordinary neutrino masses 
assures that the mirror fermion contributions to LFV processes are by far the dominant 
effects. Moreover, the absence of QCD corrections and hadronic matrix elements allows 
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in most cases to make predictions entirely within perturbation theory. Exceptions are 
the decays B ds — > £i£j that involve the weak decay constants Fs d a ■ 

The decays and transitions considered by us can be divided into two broad classes: 
those which suffer from some sensitivity to the UV completion signalled by the logarith- 
mic dependence on the cut-off (class A) and those which are free from this dependence 
(class B). We have 

Class A 

H~ — > e~e + e~ , t~ — > e~e + e~ , t~ — > fi~fi + fi~ , ji — )■ e conversion , 
— t- fi~e + e~ , r~ — t- e~/i + /i~ , r— > {i(e)n , r— fi(e)rj , r — >■ fi(e)rj' . 

Class B 

fi — > ej , t — y cy , r — > fij , 
K L ,s -> Ate , K L)S ->■ it fie , B d)S fie , B dtS re B diS r/i , 

e"/i + e" , r~ /i"e + /i~ , ~ 2) M . 

Clearly the predictions for the decays in class B are more reliable, but we believe that 
also our estimates of the rates of class A decays give at least correct orders of magnitude. 
Moreover, as pointed out in [3], the logarithmic divergence in question has a universal 
character and can simply be parameterized by a single parameter <5di v that one can in 
principle fit to the data and trade for one observable. At present this is not feasible, 
but could become realistic when more data for FCNC processes will be available. This 
reasoning assumes that <5di v encloses all effects coming from the UV completion, which is 
true if light fermions do not have a more complex relation to the fundamental fermions 
of the UV completion, that could spoil its flavour independence. 

Bearing this in mind the main messages of our paper are as follows: 

• As seen in Table 2, several rates considered by us can reach the present experimental 
upper bounds, and are very interesting in view of new experiments taking place 
in this decade. In fact, in order to suppress the /i — > e~f and — > e~e + e~ decay 
rates and the fi — e conversion rate in nuclei below the present experimental upper 
bounds, the relevant mixing matrix in the mirror lepton sector, Vne, must be rather 
hierarchical, unless the spectrum of mirror leptons is quasi-degenerate. 

• The correlations between various branching ratios analyzed in detail in Section 13 
should allow a clear distinction of the LHT model from the MSSM. While in the 
MSSM without significant Higgs contributions the dominant role in decays with 
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three leptons in the final state and in [i — e conversion in nuclei is played by the 
dipole operator, this operator is basically irrelevant in the LHT model, where Z°- 
penguin and box diagram contributions are much more important. While in the 
Higgs mediated case, the distinction of the MSSM from the LHT model is less 
pronounced, all ratios involving ^ — > £jj and in particular decays with electrons 
in the final state still offer excellent means to distinguish these two models. 

• The measurements of all rates considered in the present paper should allow the full 
determination of the matrix Vhi, provided the masses of the mirror fermions and 
of the new heavy gauge bosons will be measured at the LHC. 

• We point out that the measurements of BriK^ — > fie) and BriK^ — > 7r°/xe) will 
transparently shed some light on the complex phases present in the mirror quark 
sector. 

• The contribution of mirror leptons to (g — 2) M is negligible. This should also be 
contrasted with the MSSM with large tan (3 and not too heavy scalars, where those 
corrections could be significant, thus allowing to solve the possible discrepancy 
between SM prediction and experimental data [69]. 

• Another possibility to distinguish different NP models through LFV processes is 
given by the measurement of fi — > with polarized muons. Measuring the angular 
distribution of the outgoing electrons, one can determine the size of left- and right- 
handed contributions separately [86]. In addition, detecting also the electron spin 
would yield information on the relative phase between these two contributions [87]. 
We recall that the LHT model is peculiar in this respect as it does not involve any 
right-handed contribution. 

• It will be interesting to watch the experimental progress in LFV in the coming years 
with the hope to see some spectacular effects of mirror fermions in LFV decays that 
in the SM are basically unmeasurable. The correlations between various branching 
ratios analyzed in detail in Section 13 should be very useful in distinguishing the 
LHT model from other models, in particular the MSSM. In fact, this distinction 
should be easier than through high-energy processes at LHC, as LFV processes are 
theoretically very clean. 

The decays /i — > ej, t — > [iir and (g — 2) M have already been analyzed in the LHT 
model in [26,27]. While we agree with these papers that mirror lepton effects in fi — > e-f 
and r — > fin can be very large and are very small in (g — 2) M , we disagree at the 
quantitative level, as discussed in the text. 
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A Neutrino Masses in the LHT Model 

Within the LH model without T-parity, there have been several suggestions how to 
naturally explain the smallness of neutrino masses within the low-energy framework of 
the model [88-91]. Although differing in the details, they are all based on a coupling of 
the form 

Yi j (L i ) T $C- 1 U i + h.c, (A.l) 

where L l are the left-handed SM lepton doublets and C is the charge conjugation oper- 
ator. This term violates lepton number by AL = 2 and generates Majorana masses for 
the left-handed neutrinos of size 

m i:j = Yijv' , (A.2) 

with v' being the VEV of the scalar triplet $. 

Thus, in order to explain the observed smallness of neutrino masses, either Y^ or 
v'/v has to be of O(l0~ 11 ). While in the case of Yy, this appears to be an extremely 
fine-tuned scenario, in the case of v' some so far unknown mechanism could be at work 
that ensures the smallness of v' . Such a mechanism would also be very welcome from 
the point of view of electroweak precision observables. 

Indeed, (A.l) is a concrete example for the general mechanism found and discussed 
in [92] to generate Majorana neutrino masses for the left-handed neutrinos through their 
interaction with a triplet scalar field. 

One should however be aware of the fact that (A.l) explicitly breaks the enlarged 
[SU(2) x U(l)] 2 gauge symmetry of the LH model, while it is invariant under SU(2)l x 
U(1)y- Consequently, in [90] the interaction (A.l) has been encoded in the gauge- 
invariant expression 

Y^UfY^C^V + h.c, (A.3) 
where £ is an SU (5) symmetric tensor containing, amongst others, the scalar triplet <E>. 
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Now, rriij is given by 

rrnj = Yij (v> + , (A.4) 

so that ~ C(10 -11 ) is necessarily required in order to suppress also the second 
contribution. 

In the case of the LHT model, an interaction of the form (A.l) is forbidden by 
T-parity. However, one could T-symmetrize the interaction term (A. 3), leading to 



Y 



+ h.c, (A. 5) 



where Q = diag(l, 1,-1,1,1). In this way, a neutrino mass matrix 



n 2 



m a = Y njj ( A - 6 ) 

is generated. Note that this corresponds to only the second term in (A.4), as the first 
term is forbidden by T-parity. So again we are forced to fine-tune Yy to be O(l0~ n ), 
which is almost as unnatural as just introducing a standard Yukawa coupling to make 
the neutrinos massive. 

A different way to implement naturally small neutrino masses in Little Higgs models 
has been developed for the Simplest Little Higgs model in [93,94]. Here, three TeV- 
scale Dirac neutrinos have been introduced, with a small (~ 0.1 keV) lepton number 
violating Majorana mass term. Like that, naturally small masses for the SM neutrinos 
are generated radiatively. This idea can easily be applied to the LHT model. However, as 
already discussed in [94], the mixing of the SM neutrinos with the heavy Dirac neutrinos 
appears at 0(v/f), so that / is experimentally constrained in this framework to be at 
least ~ 3 — 4 TeV. This bound is much stronger than the one coming from electroweak 
precision constraints, / ^ 500 GeV [12], and re-introduces a significant amount of fine- 
tuning in the theory. For that reason, we do not follow this approach any further. 

Thus so far we did not find a satisfactory way to naturally explain the smallness 
of neutrino masses in the LHT model. In fact, it is easy to understand why this does 
not work. In order to keep the relevant couplings of 0(1), there should be either a 
very small scale, as v' in the LH model, or a large hierarchy of scales, as in the see-saw 
mechanism [95], present in the theory 9 . However, in the LHT model the only relevant 
scales are 

v = 246 GeV , /~lTeV, 47r/~10TeV, (A.7) 
which are all neither small enough nor widely separated. 



9 Note that also the smallness of the lepton number violating scale in the model of [94] would require 
some explanation. 
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To conclude, we find that the LHT model cannot naturally explain the observed 
smallness of neutrino masses. This, in our opinion, should however not be understood 
as a failing of the model. After all, the LHT model is an effective theory with a cutoff 
of O(10 TeV), while the generation of neutrino masses is, as in see-saw models, usually 
understood to be related to some much higher scale, which will in turn be described by 
the UV completion of the model. Consequently, in our analysis we have simply assumed 
that the mechanism for generating neutrino masses is incorporated in the (unspecified) 
UV completion, and that the details of this mechanism have only negligible impact on 
the low-energy observables studied in the present paper. 



B Relevant Functions 



D (x) 
E (x) 
D' (x) 



E' (x) 



4 , -19x 3 + 25x 2 x 2 (5x 2 - 2x - 6) , 

'9 lQgX+ 36(a: -1)3 + 18(s-l)* ^ 

2, x 2 (15 - 16x + 4x 2 ) , x(18-llx-x 2 ) 

-- lOg X H ; — log X H ; — 

3 5 6(1 -x) A B 12(1 -x) 3 

3x 3 - 2x 2 , 8x 3 + 5x 2 - 7x 

■ log X + 



2(x - ly 



12(x - l) 3 



3x 2 x 3 — 5x 2 — 2x 

logx - 



2{x- 1) 



4(x - l) 3 



MVi) = ~ 
= - 



V% log Hi 



I - yt) 2 1 - yi 

yf log Vi i 



l - yi) 2 l - yt 



F un (y t , z; W H ) = - Vi - F 5 (y h z) - 7F 6 (y u z) - 9U (y u z) , 
F di [y h z- W H ) = lyi - F 5 (y u z) - 7F 6 [y h z) + 3U {y h z) . 



F$ (yi,z) 
U (yi,z) 



yf log yt 



z 3 log z 



(1 - yt) (z - y^ (1 - z) (y t - z) 



y\ log yt 



+ 



z 2 log z 



(1 - y { ) (z - y^ (1 - z) (yi - z) 



yf^ogyi 



+ 



z 2 log z 



+ 



(y t -z)(l- Vl f (z - y t ) (1 - z) 2 (1 - Vi ) (1 - z) 



(B.10) 
(B.ll) 
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G (y h z; Z H ) 
Gi {y[,z'-A H ) 
G 2 (yu z\ rj) 



--U (y h z) 



25a 



G (y'i, z- Z H ) 



yf log Vi 



+ 



10a (1 - yi) (77 - yi) (yi - z) 
z 2 log z 



+ 



rf log 77 



(1 - z) (rj - z) (2; - y t ) (1 - 77) (yi - 77) (77 - z ) 



(B.12) 
(B.13) 



(B.14) 



Codd(i/0 

-Dodd(?/i) 

S'odd 



64 T 2 

4/3 



y^odd - 8yiR 2 (yi) + + 2yiF 2 (yi) 



D ( Vi ) - 7 -E ( yi ) - ^E (y[) 



log 



M w H 



(B.15) 
(B.16) 
(B.17) 



F( Zi , yj ;W H ) 



Ai(zi,yf, Z H ) 



{i- Zl ){i- yj )\ l h Vj 



zf log Zi 



y ] log yj 



0* -Vi)' 



1 - 2z 4 



[zi - y 3 )(l - Zi 



Ziyj 



1 - 2 Vj 



100a 



z' iZi log 4 



(1-m-y') (zi-y^l-ztf 



2/y2/i log 2/' 



( Zi - yj )(l-ytf\ 



3 



log a 



«i log Zi 



(a - 1)(1 - zj)(l - y') - y s )(l ~ zt)(l - z[) 



(zi-y^l-y^l-y'A 



^yj ' 



(B.H 



(B.19) 



(B.20) 



L 2 ( Vi ) 



1 



12(1 — yi) 4 
1 

6(1 



[-8 + 38y< - 39y 2 + 14yf - 5yf + 18y 2 logj/J , (B.21) 
-10 + 43y, - 78y 2 + 49^ 3 - Ay 4 - l&yf log ft] . (B.22) 
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